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ABSTRACT 



Numerical prediction is rapidly becoming the most accurate ap- 
proach to forecasting tropical cyclones. The numerical model in this 
study utilizes the U. S. Navy Fleet Numerical Weather Facility's (FNWF) 
so-called decomposition (SR) fields to produce geostrophic steering 
currents for tropical storms and/or hurricanes and typhoons. 

Fifteen of the 1965 hurricanes and typhoons were used to obtain 
over 200 twelve and twenty-four hour storm predictions, in one-hour 
time steps, for each of the following FNWF analyses: 1000, 500, 200, 
1000/500, 1000/200, and 500/200 -mb SR fields. 

Results indicate that the 500-mb SR fields, including a statisti- 
cal correction for latitude and/or longitude, yielded the least fore- 
cast error for Pacific typhoons with average 12 and 24-hour storm 
errors of 69 nm and 154 nm, respectively. Forecasts for Atlantic hur- 
ricanes verified best for 12-hour forecasts at 1000 mb, using a similar 
geographical modification, with an average error of 97 nm. For 24 hours 
the hurricanes were steered best at 500 mb, with modification, result- 



ing in an average forecast error of 177 nm. 
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Vgi 

Vgj 



height of constant pressure surface 
gravity 

the pressure altitude in the standard atmosphere 
the Coriolis parameter, 2 (a) sin Q 
grid distance, 381 km at 60N 
map scale factor 

del operator on a constant-pressure surface 

horizontal Laplacian operator on a constant-pressure 
surface 

component of geostrophic wind along the i axis of the 
FNWF grid 

component of geostrophic wind along the j axis of the 
FNWF grid 



SR-( ) FNWF decomposition field at ( ) mb. 

SR-( )/( ) FNWF decomposition field for layer between ( ) and 

( ) mb 

vector geostrophic wind 

A y\ 

i, j unit vectors along cartesian coordinate axes, I, J, of the 
FNWF grid 

x, y scalar coordinates coincident with I, J axes. 
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Introduction 



Accurate prediction of tropical storm tracks is of vital interest 
to the military services and the general public. This interest has 
manifested itself by the many subjective and objective forecast tech- 
niques developed in the past ten years [1]. Prominent among them is 
the statistical model produced by the Travelers Research Center [2], 
the NHC64 model developed by the National Hurricane Research Labora- 
tory [3] and the tropical steering program used by the United States 
Navy Fleet Numerical Weather Facility, (FNWF)[4J. The purpose of this 
study Is to continue the search for an optimum forecast technique by 
utilizing numerically produced tropospheric "long wave patterns" to 
generate geostrophic steering currents for the tropical storms. Both 
typhoons and hurricanes from the 1955 storm season in the western North 
Pacific Ocean and the North Atlantic Ocean were used in this study. 

The numerical fields were supplied by FNWF . 

All computations were made on the U. S. Naval Postgraduate School 
Control Data Corporation 1604 digital computer. 



The term "tropical storm is used throughout to generally indicat 
both the storm and hurricane /typhoon stages . 
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2 . 



FNWF Decomposition Fields. 



The so-called decomposition fields used in this study, and cur- 
rently used by FNWF in tropical storm forecasting, are isobaric 
height fields decomposed into additive components by applying a mathe- 
matical smoothing operator to the original height pattern for a pre- 
scribed number of times. [5], [4]. The smoothed height field, called 
the SR or residual analysis, simulates a long-wave circulation pattern 
which tends to be quasi-stationary . If the SR field is subtracted from 
the original height analysis, a small-scale (i.e., short wave) dis- 
turbance field (SD) is produced. 

Z = SR + SD (1 ) 

Tropical cyclones are examples of SD components. Figures 1 and 2 show 
examples of SR analyses at 1000 and 500 mb, respectively. 
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3 . 



Synopsis of the Present Technique of Numerical Storm Prognosis 
at FNWF. 



FNWF makes predictions of tropical storm movements upon request 
from the major Fleet Weather Centrals [4]. The Weather Centrals pro- 
vide FNWF with the storm position and FNWF, in turn, produces a numeri- 
cal forecast of the storm movement. The Weather Centrals use the 
numerical product as an aid in formulating the official forecasts 
which are distributed to the Fleet users [6], [7j. 

The numerical program which generates the steering trajectories 
is called HAT. The FNWF programs, including HAT, use a linear I, J 
grid, instead of latitude and longitude, to locate a point on earth [4] 
The grid is a 63 x 63 square and is centered on the North Pole. There 
are 3969 intersections in the grid with a spacing of 381 km at 60 N. 

An arbitrary nine-^)oint subset of the FNWF operational grid is shown in 
figure 3. 

HAT presently uses a space mean geostrophic wind as the basic 
steering parameter. The influence of the storm is removed from the 
observed flow by a particular smoothing operation [4], the result being 
the large scale geostrophic flow. The storm is assumed to be a point 
vortex and to have no interaction with the large scale steering current 

The mean geostrophic steering wind has the form 



and is determined at two levels, 1000 mb and 500 mb, at one-hour inter- 
vals up to 72 hours. At the 500mb level HAT uses hourly prognostic 
fields, smoothed in the prescribed fashion [4], to establish the mean 




( 2 ) 



wind. The low-level steering uses the most recent 1000-mb SR analysis. 
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HAT utilizes long term climatology in addition to the current data. 
As a result, three individual storm trajectories are produced. The 
first uses 500-mb data, the second uses 1000-mb SR analysis, and the 
third is an empirical combination of the first two and climatology. 



12 



4. 



Selection of Cases and Data. 



Data fields were obtained from FNWF for the period beginning 15 
August 1965 and ending 1 November 1965. All five 1965 hurricanes 
and ten typhoons occurred during this period. The storms were steered 

■N 

with the SR-1000, SR-500, SR-200 and the three corresponding thickness 
fields, 1000/200, 1000/500 and 500/200 mb. 

The U. S. Fleet Weather Central/Joint Typhoon Warning Center, 

Guam (FWC/JTWCG) and U. S. Fleet Weather Facility, Jacksonville, Florida 
(FWFJ) , are principally involved with tropical storms and each publish- 
es a post analysis in the form of annual storm reports [6 and 7]. These 
reports contain tabulations of tropical storm positions for one calendar 
year. One of the many entries is the Best Track position. The Best 
Track positions are determined during post-season research and can be 
generally defined as a smooth line running through the fixes. The fixes 
are determined by aircraft reports, radar, and weather satellite pictures. 
The displacement of the Best Track from the fixes varies with the method 
of determining the fixes. In most cases the displacement lies within 
10 to 20 nm [6] . 
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5. 



Procedure . 



This author wrote a numerical tropical storm steering program 
using FNWF's HAT program as the basic model. The principal changes 
made to the existent basic program allowed rapid cycling of many differ- 
ent storms on one computer run and the elimination of all empirical 
operations. In addition, only one SR field was used on a given computer 
operation. 

The Best Track positions of the tropical storms were accepted as 

2 

the initial and verifying positions. The procedure used an SR analysis 
field to obtain the storm's steering current for a twelve-hour period. 

As an example, for the 00Z upper air observation^ the 18Z Best Track 
position of the previous day was chosen as the initial storm position. 
The storm was then moved forward using 12 one-hour time steps. The 
error was determined by verification at the time of the 06Z Best Track 
position. The procedure was initialized with either the 00Z or the 12Z 
analysis and was iterated for all verifiable positions. The starting 
time and the number of iterations were dependent entirely upon the in- 
dividual storm. The model was later modified for 24-hour prediction by 
increasing the number of one-hour time steps to 24, but otherwise pro- 
ceded as outlined above. 

Geostrophic steering components were determined using the finite 
difference form of the wind equation 




The sine function in the coriolis parameter, £ , was modified at low 
2 

Best Track positions in the Atlantic included some cases of 
tropical cyclones in the depression or extratropical stage. 
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latitudes to approach that at 7 . 2 degrees in order to avoid division 
by a number approaching zero. See figure 4. 

Upon reading the initial position of the storm, the computer pro- 
gram translated latitude and longitude into an I,J position in the grid. 
From the one I,J position four surrounding positions were determined. 
These four points were as depicted in figure 3; 1+1, J; I-1,J; I,J+1; 
and I,J-1. 

An interpolation was made in the component geostrophic wind fields 
( \/gi, \Vgj) to determine their respective values at each of the four 
grid points. The one-hour movement for each of the four points was 
simply the vector combination for the Wgi and >ygj components in 
knots. Each of the four points was then moved forward one hour and the 
centrum of the new positions was defined as the new "zero plus one hour" 
I,J position of the storm. The I,J values were reconverted back to 
latitude and longitude after the 12- or 24-hour prognostications. 

Correlation and regression coefficients between true positions and 
prognostic positions were obtained in order to correct bias in the fore- 
casts. Using the statistical corrections thus obtained, the storm pre- 
dictions were then re-run for selected fields. 
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6 . 



Results for Twelve-Hour Forecasts. 



Table 1 lists the 12-hour forecast errors for each storm and the 
overall averages. The method giving the least error for a given storm 
is noted in the table by the asterisked values. The error for the SR- 
200 is highest with an overall average of 186 nm for a 12-hour movement. 

In addition, the thickness fields, using the SR-200, were also poor, 
with a 143 nm error for the SR 200/1000 and a 169 nm error for the SR 
200/500. 

The SR-1000 and the SR 500-1000 were nearly equal with overall 
error values of 108 nm and 105 nm, respectively. The SR-500 was over- 
all the best with an average error of 89 nm. While the SR-500 did not 
always produce the least error it was, nevertheless, uniformly good for 
all storms. 

After the hurricane and typhoon data were separated, it was interest- 
ing to note that the SR-1000 was best for hurricanes (Elena was an ex- 
ception) , with an average error of 101 nm. Comparing the performance 
of SR 500 and 1000 mb, SR-500 was best for typhoons (Lucy an exception), 
with an average error of 77 nm. 

A plot was made of all storms but only two are included in this 
test. (See figures 5,6,7 and 8). It was noted that the predicted posi- 
tions, especially north of 25-30 degrees, fell to the south of the 
verifying position in nearly every case. The shortcoming of under- 
forecasting poleward storm movements is consistent with other numerical 
techniques [1]. However, an empirical correction was readily derived 
to improve forecast accuracy, as follows. Regression equations re- 
lating predicted latitude to true latitude and predicted longitude to 
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true longitude were derived using both the SR-1000 and SR-500 data. 
Inasmuch as there were 91 hurricane and 115 typhoon forecasts, the 
storms were stratified by ocean as well as grouped together. The 
pertinent regression equations are listed in table 3. 

No single empirical correction was found to uniformly reduce all 
hurricane or typhoon prognostication errors. With reference to the SR 
100-mb forecast verification (see table 1), a latitude correction ap- 
plied to typhoons responded best with an overall average error reduction 
from 111 nm to 99 nm, while the longitude modification reduced the aver- 
age error by only 2 nm to 109 nm. For the hurricanes, the latitude 
correction lowered the error by 1 nm to 100 nm while the longitude 

equation decreased the error to 97 nm. A combined latitude and longi- 

tude correction to hurricanes still produced an error of 97 nm. 

The SR-500 forecasts were modified by a linear regression correc- 
tion only when the prognostic position was north and east of the init- 
ial position. Twenty-five twelve-hour typhoon forecasts were affected 
by this change. The average forecast error for all typhoons was there- 
by reduced from 77 nm to 69 nm. Typhoons Mary and Rose moved north and 

west only and were not changed. The modification, in the case of 

typhoons Della and Lucy netted an increase in average error of 23 nm 
and 7 nm, respectively. The regression technique, using the SR 500 mb 
data, raised the average error for hurricane forecasts from 101 to 121 
nm. 

The SR steering method shows improvement over the official fore- 
casts in three out of five hurricanes as may be seen in the last column 
of Table 1. The average official forecast and modified SR 1000 errors 
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are comparable, 94 and 97 nm, respectively. Twelve*.hour official fore- 
cast errors for typhoons were not available for comparison. 

The major reason that the numerical steering model, using the SR- 
1000 data did not produce a lower average, especially in the Atlantic, 
was the model* s failure to predict the high latitude movements of 
hurricane Elena. Elena became a hurricane at about 23° north and con- 
tinued as a hurricane to about 45° north. Typhoon Carmen and hurricane 
Carol also share similar problems when poleward of 25-30° latitude. 

The reader will note the similarity in track orientation, if not 
position, in figures 5,6,7, and 8. When Carol made a loop in mid- 
Atlantic both the SR 1000 and the SR 500 made similar loops (figures 5, 
6) even though displaced to the southeast. Moreover, both the SR 1000 
and the SR 500 predicted the completion of the looping evolution. 

Such behavior is not unusual in the cases studied. 
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7. Results for Twenty -Four Hour Forecasts. 

Table 2 lists the storms and composite 24-hour forecast errors for 
both official forecasts and the SR approach. Only the SR 500 was given 
extensive study for the 24-hour forecast because of the limited amount 
of time available. The SR 500 was chosen because of its good perfor- 
mance on the 12-hour predictions and because it out performed the SR 
1000 and the SR 200 on the two 24-hour test cases, Carol and Carmen. 

The SR 500 predictions, unmodified by latitude and/or longitude 
corrections, were slightly higher than the official forecasts for both 
hurricanes and typhoons. 

The 24-hour forecasts for SR 500 were modified for both typhoons 
and hurricanes as a single group in a manner similar to the 12-hour 
SR 500 predictions. However the modification statistic was derived 
only for those storms forecasted to move towards the northeast. 

All hurricane prediction errors decreased with the modification 
except for Anna. The average was reduced from 198 nm to 177 nm. The 
average storm error, using official forecasts, was 187 nm for hurricanes 
[7]. 

The prediction error for typhoons by storm was reduced from 168 nm 
to 154 nm with the modification. Two typhoons, Lucy and Olive, showed 
increased average forecast errors of 20 nm and 12 nm, respectively, when 
modified. Again, as in the 12-hour predictions, typhoons Mary and Rose 
did not move toward the quadrant east of north and, therefore, were not 
modified. The average official forecast error for typhoons was 153 nm 
[ 6 ]. 

The storm tracks (figures 4 and 6) for the 24-hour SR 500 follows 
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the Best Track quite well. When Carol began her mid-Atlantic loop, 
the 24-hour prediction shows a similar but nondescript trajectory. 

The completion of the looping action for Carol on 25 September 1965 
was forecasted very well with both the actual and forecast positions 
indicating definite termination of the recurvature path. 

The 24-hour forecasts also under predicted northward movement for 
storm positions poleward of 25° north, similar to the 12-hour forecasts. 
The magnitude of the error was greater than for the 12-hour predictions 
as can be seen in the regression equations. (See table 3.) 
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8 . 



Errors Inherent in Computations. 



Both the initial and the forecast positions of each storm were 
printed by the computer program. An unexpected result came after check- 
ing the initial position for accuracy. Of the 122 positions given de- 
tailed checks, 71 were correct, 46 were in error by 0.1 degree of 
latitude or longitude, and 5 were in error by a 0.1 degree on both 
latitude and longitude. No pattern or position bias was found for the 
errors. 

Both the data cards and the basic computer program were found to 
be correct. The explanation for the loss of accuracy is in the conver- 
sion of latitude and longitude to I , J values and the reconversion from 
I, J back to latitude and longitude. An insufficient number of Binary 
Digits are allotted in the sub-routines used for the fractional por- 
tions of both I, J values and latitude and longitude in the conversion 
seouences. For example only four Binary Digits are used in the frac- 
tional portion of the conversion of I, J to latitude and longitude. A 
minimum of seven Binary Digits would be required to insure one decimal 
place accuracy. 
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9. 



Conclusions . 



a. It is concluded that the 12- and 24-hour movement of 
tropical storms can be predicted using the FNWF de- 
composition analyses (SR fields) with accuracy com- 
petitive to that of operational techniques now in exist- 
ence . 

b. Secondly, it appears that typhoon movement can be pre- 
dicted more accurately than hurricane movement. 

c. Further, the observed 12-hour hurricane motion was found 
to correlate better with the SR-1000 while the 24-hour 
movements correlated better with the SR-500. 

d. In addition, the SR-500, modified for latitude bias, 
improves the typhoon forecasts for both 12 and 24-hours. 

Several facts temper the success noted above. One, the use of 
Best Track position here gives the official forecast a relative handi- 
cap, and two, the number and time of the experimental forecasts differ 
from the official forecast sample. Thirdly, it is to be noted that the 
SR field employed represents the analysis at the midpoint of the 12- 
hour forecast interval while the 24-hour forecast used the analysis 
one-quarter of the way from initial to verification time. Other more 
exacting tests may now be employed as the subject research suggests 
success in the realm of a feasibility study. 
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10. Recommendations . 



1. The numerical geost rophic steering should be tried in a 700- 
mb SR field and the results compared to the SR-1000 and SR-500 results. 
The 7 00 -mb level is used in several current prediction models. In 
addition, some other intermediate level between the 1000 -mb and 500- 
mb levels may well be superior. 

2. A vertically integrated wind over the entire atnosohere ought 
to bo given extensive testing as a steering current. 

3. The prognostic SR fields for 24 and 48 hours may be the media 
for improving tropical storm extended forecasts. 

4. The number of binary digits should be increased in the frac- 
tional portions of the conversions and re-conversions of latitude and 
longitude to I, J, in order to improve accuracy in tropical storm fore- 
casts . 

5. Publication of 12-hour official forecast errors is most 
desirable, since this time oeriod should be considered fully prior to 
attempting extended predictions. 
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Overall Average: 164 178 
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Hurricanes and typhoons moving with component toward the north and east: SR 500, latitude, 24-hour Coreca 
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